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Abstract 

We discuss whether some perturbed Friedmann-Robertson- Walker (FRW) universes could be 

Q-i. 

<^ . creatable, i. e., could have vanishing energy, linear momentum and angular momentum, as it could 

^O | be expectable if the Universe arose as a quantum fluctuation. On account of previous results, the 

q ' background is assumed to be either closed (with very small curvature) or flat. In the first case, fully 

i ■ arbitrary linear perturbations are considered; whereas in the flat case, it is assumed the existence 

of: (i) inflationary scalar perturbations, that is to say, Gaussian adiabatic scalar perturbations 

having an spectrum close to the Harrison-Zel'dovich one, and (ii) arbitrary tensor perturbations. 

We conclude that, any closed perturbed universe is creatable, and also that, irrespective of the 

OO ■ 

spectrum and properties of the inflationary gravitational waves, perturbed flat FRW universes with 

standard inflation are not creatable. Some considerations on pre-inflationary scalar perturbations 

OO 
O 



are also presented. The creatable character of perturbed FRW universes is studied, for the first 
time, in this paper. 
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I. INTRODUCTION 



In |l| (see also the addenda in [2J), for a wide set of non asymptotic Minkowskian space- 
times, it was uniquely determined when one of these space-times can be said to have van- 
ishing energy, vanishing linear 3-momentum and vanishing angular 4-momentum, that is, 



vanishing energy and momenta. These vanishing values cou 



d be expected in the case of a 



universe which rose from a quantum vacuum fluctuation |3j [4|. In [l|, universes of this kind, 
with vanishing energy and momenta, were called 'creatable universes'. 

It is very well known that, whatever the energy-momentum complex may be, the definition 
of energy and momenta of the universe is strongly dependent on the coordinate system. One 



11 



must then stress the uniqueness reached (l| in the characterization of the family of creatable 
universes, within the above wide set. This uniqueness has been reached by using some 
physical criteria to select the appropriate coordinate systems. See [5| as an example of a 
more mathematical criterion, which leads the author to use conformally flat coordinates in 
the space-time, when it is possible. 

In order to look for convenient physical criteria to select the appropriate coordinates, 
one should assume that the proper energy and momenta of any space-time representing 
the universe is conserved in time. Taking into account that the Universe is supposed to 
embrace everything, fluxes of energy and momenta going out of such an entity are not 
possible and, consequently, the above assumption about time conservation seems actually 
appropriate. Such an assumption only has a clear physical meaning if it is referred to a 
physical (proper) and universal (synchronous) time. In other words, the conserved energy 
and momenta of a universe should be calculated by using Gauss coordinates, which involve 
the required time coordinate. Energy and momenta would be then conserved in terms of 
Gauss time coordinates. Because of this conservation, no global Gauss coordinates are 
necessary in practise to calculate energy and momenta. Actually, Gauss coordinates defined 
in the elementary vicinity of a generic space-like 3-surface suffice. 

Given a space-time, there are plenty of different Gauss coordinate systems. Thus, we will 
need to make sure that the characterization of any perturbed FRW universe as creatable 
or non creatable is independent on the Gauss coordinate used. We will consider various 
of these perturbed universes and, then, in order to see if they are creatable, we will apply 
different strategies, inspired by [if, which are adapted to each particular case. Nevertheless, 



any of these strategies will obey the following protocol: (i) take any space-like 3-surface, S3, 
and build the corresponding Gauss coordinates in its elementary vicinity, (ii) look for new 
Gauss coordinates leading to an 'instantaneous' 3-space metric, (11% = dl^t = t ), which 
explicitly exhibits its conformal character on the boundary, E 2 , of £3 (these coordinates 
always exist; see |l|) an d> ( m ) calculate the energy and momenta of the universe in the 
resulting coordinate system. Despite the still remaining freedom in the election of these 
coordinates, it can be seen that the calculated values of energy and momenta are unique 
for the different cases we consider in the present paper. In all these cases, if these energy 
and momenta vanish, the universe is creatable, since it can be directly seen that they vanish 
irrespective of both the selected space-like 3-surface, £3, and the chosen time, to (energy 
and momenta conservation). This will be verified at the end of each particular case. 

It has sometimes been argued (see, e. g., J5|) that, in a space-time which is not asymp- 
totically flat, the global energy and momenta would have no physical meaning, since these 
quantities could never be measured; however, it is also claimed that the energy and mo- 
menta of any part of this space-time have physical meaning. This double claim is not fully 
consistent; in fact, if energy and momenta have physical meaning for any part of the space- 
time, the corresponding global quantities should be interpreted as the limit of the physically 
significant energy and momenta of the parts, as they grow to fill the entire space-time. 

By using the Weinberg energy-momentum complex [10J and the above protocol, it was 
proved that the closed and flat Friedmann- Robert son- Walker (FRW) universes (K = 0, +1) 
are creatable, whereas the open version (K = —1) is not [l|]. This is in agreement with most 
of the papers on the subject, but not with all of them (see again |5|). Similar conclusions 
were obtained by using very different methods. For example, in papers 6j and [7[, it was 
proved that the tunneling amplitude for creation from 'nothing' is finite in the case of closed 
cosmological backgrounds and for a flat De Sitter unperturbed universe. These two papers 
concluded that the tunneling amplitude vanishes in the open case. These conclusions were 
obtained in the framework of particular FRW models for tunneling. Realistic universes with 
perturbations were not considered at all. The main goal of the present paper is the study of 
perturbed FRW universes without any calculation of transition amplitudes. Instead of these 
calculations, we use the above protocol to look for creatable universes (those with vanishing 
energy and momenta) without modelling any quantum tunneling. 

The observations seem to indicate that we live in a perturbed FRW universe and, con- 



sequently, if the above protocol makes sense, at least one of the perturbed FRW universes 
should be creatable. Our results confirm this expectation, since we have proved that any 
perturbed closed universe is creatable. The perturbed K — — 1 universe does not need to be 
considered here because the corresponding background is not creatable. The perturbed flat 
case is studied in detail along the paper. It is not creatable, under very general conditions 
strongly supported by current observations. 

This paper is organized as follows: In Sec. [Til we consider absolutely arbitrary linear 
perturbations in the case K — +1. The case K = is studied in Sec. IIHI The total 
energy due to the scalar perturbations arising in standard inflationary models is calculated 
in Sec. lHI Al (it is infinite). Sec. IIII Bl contains the calculation of the same quantity in the case 
of fully arbitrary gravitational waves (it vanishes). Finally, in Sec. IIV[ we summarize our 
main conclusions and present a general discussion, including comparisons between our results 
and methods and those of previous papers dealing with quantum creation from 'nothing' 
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Some of these results have been briefly presented, with no calculations, in the meeting 



ERE-2007 [11] 



Let us finish this section with some words about notation. Units are chosen in such a 
way that the speed of light is c = 1. The gravitational constant, the scale factor, and the 
index of the 3-space curvature are denoted G, a, and K, respectively. Symbols t and r stand 
for the Gauss and the conformal times (adr = dt) and, finally, the unit vector k and the 
modulus k define a generic vector k = kh in momentum space. 

II. THE CASE OF A PERTURBED CLOSED FRW UNIVERSE 

The line element of the closed FRW space-time can always be written in the form: 

ds 2 = -dt 2 + dl 2 , dl 2 = — — — ^S ij dx i dx j , r 2 = S ij x i x j (1) 

(1 + f r 2 ) 

with K = +1. Here, x l and t are global Gauss coordinates. Moreover, the 3-space metric 
exhibits a conformally flat form everywhere. 

Now, let us consider the case of a perturbed closed FRW universe with scalar and tensor 



perturbations [12]. Then, in the synchronous gauge (see the Appendix), conditions goo = — 1, 



<7io = are satisfied. Hence, in this gauge, Gauss coordinates are used. In terms of them, 



the 3-space metric, dl , reads as follows: 

dl 2 = — a ® 2 (5ij + hij)dx l dx\ (2) 

(1 + \r 2 ) 

where the hij(t,x l ) functions are such that hij << 1. 

Let it be a particular space-like 3-surface, t = to, and particularize the above 3-space 
metric on this 3-surface, that is, consider dl(t = t ) 2 = dl^. This instantaneous 3-space 
metric is a conformally flat 3-metric on the boundary 2-surface, S 2 , of t — to (l|. Then, 
among the different Gauss coordinate systems, one always can select some ones such that dl^ 
on the boundary £2, say dlo\j: 2 , shows explicitly its conformally flat character. According to 
the protocol displayed in Sec. [H which can be used here to calculate the energy and momenta 
of perturbed closed universes, we can pick up any one of these last coordinate systems to 
compute the corresponding values of the energy and momenta of our perturbed universe. 
We will subsequently see that these values are unique. 



Then, according to 10[, the energy, P°, the 3-momentum, P l , and the 4-angular momen- 



tum, (Ji k , J 0t ), of the universe are: 

P° = ^J(d jgij -d ig )dE 2i , (3) 

p% = T^oJ^-gij)^^, (4) 

J 3k = 7^-77 / ( x k9ij - Xjgki)dE2i, (5) 

l07TG J 

J 01 = P l t - ^ J [(d k9kj - d 3 g)x t + gS tj - g tJ ]dZ 2j , (6) 

where g\j = dtgij, g = dijgij, and where dY>2i stands for the integration element on £2. In 
the present case, the 2-surface £ 2 is r = 00. 

Now, it is obvious that the energy and momenta of our perturbed closed FRW universe 
must vanish; in fact, according to Eq. (J2J), all the integrands in the above 2-surface integrals 
go at least like 1/r 4 as r tends to 00. (Notice that h^ cannot grow indefinitely with r, 
otherwise we would not have h^ << 1 everywhere). Of course, this asymptotic behavior, in 
terms of 1/r, is valid for all Gauss coordinates which preserve the perturbed character of the 
metric Eq. ([2j), and so, whatever be the time parameter to, it is valid for the particular Gauss 
coordinate system where the 3-metric, dl^, shows explicitly its conformal flat character on 



E2. Evidently, such a behavior of the integrands implies straightforward that all the integrals 
in Eqs. (JSJ)-© vanish irrespective of the above Gauss coordinate used. Consequently, we 
can state that, whatever the properties of the perturbations may be, perturbed closed FRW 
universes are creatable. 

III. THE CASE OF A PERTURBED FLAT FRW UNIVERSE 

In the case of a flat background, we cannot invoke the strong decaying of the 3-space 
metric, as r tends to oo, to conclude that energy and momenta vanish; namely, arguments 
similar to those of the last section do not apply. Furthermore, one could erroneously think 
that, in the flat case, the vanishing of the energy would come straightforward from the 
cosmological principle; namely, from the assumed statistical homogeneity when averaging 
on large enough volumes. Actually, as we will see in detail along this section, the values of 
the energy and momenta of the universe depend on both the statistical character and the 
spectra of the perturbations. 

A. Scalar perturbations 



Let us first consider the case of scalar perturbations 12J in a flat FRW universe. In any 
synchronous gauge (Gauss coordinates), the perturbed 3-space metric is gr^ = a 2 (r)(6ij + hij) 



and, according to [13j, the metric perturbation, hy, can be expanded in scalar harmonics 
(plane waves) as follows: 

Mx, t) = ! d'ke^hijik, r)= ( d 3 fce ik - x [k l k J 7i(k, r) + (k*tf - ^ ij )6r/(k, r)] , (7) 

where two functions, /i(k, r) and r/(k, r), defined in momentum space, have been introduced. 
Notice that, in order to have a real value for hij, functions h and r\ must satisfy the conditions 
h(-k) = -h*(k) and rj(-k) = -rf(k). 

According to the above protocol to calculate energy and momenta, we must use a new 
Gauss coordinate system, in which, the transformed components of the instantaneous 3-space 
metric, g^At = to), show explicitly its conformal flat character on the boundary r = oo. In 
the present case, we can do more than this since our instantaneous 3-space metric is a 
conformally flat metric everywhere on the space-like 3-surface t = to ( n ot only on the 2- 
surface r = oo). In order to prove this statement, let us work in the k-space. Given a 
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k-value, let us consider the infinitesimal coordinate transformation 

x 4 = x'* + e*(x,k), (8) 

which leads to the following equation: 

h'^y^ = M k K k '* + d i e J + d J e i ■ ( 9 ) 

If we then choose e l (x, k) = /(k)e lkx k J , where /(k) is defined by the relation 

(h + Q V + 2a 2 k 2 f) t=to =0, (10) 

equation Q leads to the following new Fourier component /^(k) for t = t 

^•(k) t=t0 ee ^.(k) = -2r /0 (k)^, (11) 

where r/o(k) = r/(k, t = to)- Equation fill I) implies that the instantaneous metric <7y(£ = to) 
is conformally flat everywhere on the 3-space t = to, as we wanted to prove. 

To calculate the energy, we must insert the new components of the 3-space metric g[At = 
to) = CLol^ij + Kjit = t )] in the integral of Eq. (J3j), which must be performed on the 2- 
surface r = oo. This integral is to be calculated in the new coordinates x h ; however, in 
practice, the old coordinates, x l , can be used in first order calculations; in fact, since the 
integrand in Eq. (J3]) trivially vanishes for the background metric (a flat FRW universe), 
this integrand is a first order quantity in the perturbed universe and, consequently, the 
integral can be evaluated irrespective of the coordinate system (the differences between the 
integrands in the old and new coordinates are second order quantities to be neglected in any 
linear approach). In short, to first order, metric perturbations and its derivatives can be 
calculated in terms of the old coordinates, and the old 2-surface element, dT, 2 i, can be used 
instead of the new one. Thus, according to Eq. ffTTT) . the total energy P° can be written as 
follows: 

P° = lim ^- [ dSlmdi ! dk 3 ^-*, (12) 

r-*oo 87TG J J 

where dQ is the solid angle element in spherical coordinates, i. e., dQ = sin6d8d<j), and 
where n^ = x % jr. 

It is noticeable that, according to the above expression of P°, the energy of the universe 
does not depend on the function h in Eq. (JTj). It only depends on the function rj , that is, on 
the traceless part of the perturbed metric ([7|). 



We assume statistical isotropy, which characterizes cosmological processes as, e. g., stan- 
dard inflation. Since there are no privileged directions, the power spectrum of r/(k), namely, 
the function P v (k) = (|?7(k)| 2 ) is independent on k. It only depends on k. Three /c-intervals 
can be distinguished (see the Appendix). The interval (0, k m i n ) involving pre-inflationary 
perturbations with some unknown power spectrum P r?1 (k), the (k m i n , k max ) interval with 
inflationary perturbations evolving outside the horizon, whose spectrum is P^(k), and the 
the interval (k max , oo) with inflationary perturbations which have reentered the horizon; in 
this case, the spectrum is denoted P v3 (k). 

Inflationary perturbations are all Gaussian in the evolution period under consideration 
(before entering in nonlinear processes, see the Appendix); hence, if pre-inflationary per- 
turbations are also assumed to be Gaussian (the non Gaussian case is discussed below), 
functions 77 and h, as well as the Fourier transform of any physical quantity (density con- 
trasts, pressure, and so on) are complex functions with random phases [14J for < k < oo 
and, consequently, the 77 values necessary to calculate P° -from Eq. ([3])- can be written in 
the form 



1 



1/2 



v (k) = -^=[P v (k)\ [i/i(k)+nfe(k)j, (13) 

where 7/1 (k) -M?/2(k) are the phases, which must be assigned taking into account the relation 
77 (— k) = — 77* (k) and the relations: 

( yi (k)> = (ife(k)) = , (14) 

(l/i(k)ife(lO) = 0, (15) 

(y n (k)y n (k')) = 5 3 (k-k f ), (16) 

where index n run from 1 to 2, vectors k and k' are arbitrary momenta, and the averages are 
to be performed in the set of the universe realizations. In the last equation, 5 3 is the three 
dimensional Dirac ^-distribution. Equations ( IT^1) -( IT6"1) ensure (central limit theorem) that 
the distribution of P° values corresponding to all the possible realizations of the universe is 
Gaussian. 

By using the relation 



/ 



rr-, ikx 47ri/sin£;r \ r 

dilriie = - — — cos kr k, 



tr V kr < 17) 



which can be easily obtained, plus Eqs. (TT2I) and (TT3l) . one easily finds: 

P° = -(al/2G) lim r V / d 3 k \P m {k)] (^— - cosArWk) , (18) 

where the index j number the three intervals defined above. The volume integrals in mo- 
mentum space are extended to the regions k < k min , k min < k < k max , and k > k max , in 
cases j=l, j=2, and j=3, respectively. In order to derive the last equation, we have used the 
relation r/(— k) = — r/*(k) and, consequently, the energy P° appears to be real valued as it 
must be. It is also important that, as a result of the same relation, quantity P° only depends 
on the random variable yi(k); whereas y^ (k) becomes fully irrelevant in our calculations. 

Taking into account Eq. (TlBl and the fact that (7/1 (k)) vanishes, one easily finds the 
relation (P°) = 0. Of course, the average is performed for all the possible realizations of the 
universe, which correspond to distinct y\ (k) values but to the same power spectrum. Since 
the phases have been chosen in such a way that the P° values are normally distributed, we 
must now calculate the variance, ((P ) 2 ), which fully characterize this normal (Gaussian) 
distribution with vanishing mean. From Eq. ( Tl8|) . one easily finds a formula for (P ) 2 which is 
the addition of six different terms. Three of them involve products of integrals corresponding 
to distinct j- values (hereafter crossed terms). Since two k- vectors corresponding to distinct j 
values cannot coincide, equation ffT6|) implies that, after averaging, the crossed terms vanish. 
From the remaining three terms plus Eq. (1161) one easily obtains the following variance: 

((pO) 2 ) = (a*/4G 2 )limr 2 J2 I ^^(^(^j^ " cosAx) 2 . (19) 

Each of the three terms of this last formula is positive (it is the square of one of the three 

integrals j = 1, 2, 3) and, consequently, if one of them diverges, the variance ((P ) 2 ) diverges. 

In the case of inflationary adiabatic scalar perturbations which evolve outside the horizon 

in the radiation dominated era {k min < k < k max ), the solution of the linearized Einstein 



field equations can be found in Ma and Bertschinger 13] (see the Appendix for more details). 
In the synchronous gauge, function rj evolves as follows: 

r ] = C + ^ + aCk 2 r 2 - (20) 

in this equation, a is a pure number (see the Appendix), whereas C and D are dimensionless 
integration constants. That is to say, C and D cannot depend on time, but they may depend 



on k. Notice that here, as well as in [13|, the components of x and k are dimensionless, since 
they are comoving coordinates and their associated momenta. We see that t] is the addition 
of three different modes: one of them grows like r 2 , another one is constant, and the third 
one decreases like r _1 . 

Let us now study the behavior of the integrals involved in Eq. (T031) as r tends to oo. We 
begin with one of the contributions to the j '• = 2 integral (inflationary Gaussian perturbations 
with super-horizon sizes). It is the contribution due to the constant mode C(k) in Eq. (12131) . 
This contribution is hereafter denoted ((P ) 2 ). The spectrum P c {k) of the constant mode 
is defined by the relation 

P c {k)] 1/2 = Ak n , (21) 

where A is a normalization constant and the spectral index is n = —(3/2) — j3. According 
to Eq. (1A.10I) in the Appendix, constant /3 — (1 — n s )/2 is positive and small as compared 
to unity. From the spectrum ()2ip . one easily finds the following variance: 

((p0)2) = (vra^VG 2 ) Urn r 2+2 ^ f ^ dy(^- - cosy) 2 ) /y l+2/3 , (22) 

where y = kr. For any k, the new variable y tends to infinity as r does. On account of this 
fact, it is easily proved that the integral in Eq. ( 122]) goes just like the power r~ 2/3 as r tends 
to infinity. Hence, ((P ) 2 ) diverges as r 2 . 

Let us now consider the growing and decaying terms appearing in Eq. ( 12(31) . It is obvious 
that the contributions of these time dependent terms to the energy cannot compensate the 
infinite value of ((P ) 2 ) corresponding to the constant mode, which implies that the total 
contribution to ((P ) 2 ) due to the j '• = 2 k- interval cannot become finite. More precisely, 
from Eq. ( 1201) we find six distinct terms contributing to the j = 2 integral. The term 
used in the previous calculation does not depend on time, whereas the remaining ones are 
time dependent. By this reason, the variance ((P ) 2 ) is proportional to a 4 for the time 
independent mode, whereas it exhibits other time dependences in the remaining cases. In 
all, compensation of the resulting terms (with distinct time evolutions) to give a conserved 
finite total energy is not possible, which ensures that the entire contribution j = 2 to (P ) 2 ) 
is infinite. Moreover, this positive contribution cannot be compensated by those of the j = 1 
and j '• = 3 cases to give, finally, a finite global value for (P ) 2 ), the reason being that, as 
explained above, the j = 1 and j '• = 3 contributions are also positive. All in all, the following 
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relation holds: 

<(P°) 2 > = oo (23) 

and, consequently, we have a Gaussian statistical distribution of P° values with zero mean 
and infinite variance. Since the Gaussian probability density is 

P(P°) = l P -(P -(P )) 2 /nP°Y 2 ) (24) 

^((P ) 2 ) 

it is evident that, in our case, the probability of any particular finite value of P° vanishes 
and, as a consequence, we can say that, in the flat perturbed FRW universe under consider- 
ation, the contribution of the inflationary scalar perturbations to the total energy is infinite. 
Nevertheless, in the standard inflationary paradigm, there are also tensor perturbations, 
whose contribution to the total energy of a flat perturbed universe is calculated in next 
section. 

Let us now discuss the case of non-Gaussian pre-inflationary perturbations. These pertur- 
bations are assumed to be generated in some process fully independent on inflation and, as 
explained in the Appendix, they are assumed to be significant only in the interval (0, /c m j„); 
hence, Eqs. (Tl6l) are satisfied except for pairs of vectors k and k' whose modulus are both 
inside the interval (0, k m i n ). Therefore, there are no crossed terms in the development of 
((P ) 2 ) (see above). As in Eq. (119J1 . the variance ((P ) 2 ) is the addition of three positive 
terms. Terms j = 1 and j '• = 2 have the same form as in the Gaussian case, whereas the 
term j = 3 would be different. In this situation, the term j = 2 diverges (same discussion as 
above) and, consequently, Eq. (|23|) holds (which strongly suggests a non creatable universe). 
However, the distribution of P° values is not Gaussian in this case, and its probability 
density should be calculated for each particular non Gaussian model. Since Eq. (12^1) does 
not apply, the meaning of an infinite variance is now less clear. By this reason, in order 
to properly prove that flat universes with standard inflation are not creatable, we prefer a 
general argument, which proves that, if ((P ) 2 ) diverges in the interval (k min ,oo), namely, 
if the total energy of the inflationary Gaussian scalar perturbations is infinite, the energy 
of the universe cannot vanish whatever the properties of the non-Gaussian pre-inflationary 
perturbations may be; namely, the universe is not creatable. This is trivially proved taking 
into account that, in any admissible universe, the total energy must be a conserved quantity 
(see Sec. [T|) which vanishes in the creatable case. In fact, a vanishing energy after inflation is 
only possible if the pre-inflationary energy is infinite and it exactly compensates the infinite 
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energy associate to the scalar inflationary perturbations. However, such a universe would 
have an infinite energy before inflation and a vanishing one after this process, which is not 
compatible with the required energy conservation. 

Now, before ending this Section, we raise some comments about the consistence of the 
main result: the infinite value of the energy we have found. 

First of all, although the infinite energy has been formally obtained for a given value of 
to and a certain S3, it is obvious, by following the implementation of the protocol we have 
used, that this to value and the choice of E 3 are both arbitrary. 

Finally, we could perform a conformal coordinate transformation and still retain the ex- 
plicit conformal flat form of the 3-space metric on E 3 (remember that, according to Eq. (1111) . 
in the present case, the 3-space metric is a conformally flat one all over E 3 ). But, trivially, 
the above infinite value of the energy does not depend on the conformal transformations 
which are pertinent here: those belonging to the translation, rotation, and dilatation sub- 
groups of the conformal group. For all these coordinate transformations the energy remains 
infinite. 

B. Tensor perturbations 

In this section we are concerned with tensor perturbations evolving in a flat background. 
In such a case, the instantaneous 3-space line element is: gij = Oq(5jj + h i -). In Fourier 
space, we can write 

hl(k) = H(k,r )e tl (k), (25) 

where the quantities e^ satisfy the conditions given in the Appendix. 

Let us calculate the energy, P°, of these tensor perturbations. According to our protocol, 
calculations must be performed in a new coordinate system, in which the instantaneous 
3-dimensional metric explicitly exhibits its conformally flat character on the boundary 2- 
surface r = 00. For each value of k, let us consider a coordinate transformation of the form 
(jEJ)- Functions e*(x, k) must be chosen in such a way that, in the new coordinates x H , the 
line element on the surface r = R, where R is an arbitrary constant (at the end of our 
calculations, this constant will tend to 00), has the form 



h y (X, 



= /(n,k)e^ k - n %, (26) 

r=R 
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where h j,-(x) = h i -(k)e* k-x is the metric perturbation in the space-like hypersurface t = r 
corresponding to the fixed mode k, and where function /(n, k) is the conformal factor of 
the metric on r = R. As it has been said above, there always exists a family of coordinate 
systems in which Eq. ( )26|) is satisfied [l|. The energy and momenta can be calculated in any 
of these coordinates. With the appropiate e; functions, from Eq. (jSJ), one easily finds the 
relation 

/i%(x)=/ l y(x) + 9 i e i + a 7 -e i , (27) 

which is equivalent to Eq. (JHJ). 

The energy, P°, is calculated by using Eq. (J3J) and the metric perturbation components 
h jj(x'); nevertheless, in the linear approach we are using (see Sec. IIHAI for details), the 
following approximations can be performed: (i) write these components in terms of the old 
coordinates x l , (ii) perform the derivatives with respect to x\ and (iii) use the old 2-surface 
element, dH l 2 , instead of dH'\. Hence, we can write 

P° = (a 2 /16irG) f d j [/^.(x) - 5^ T '(x)] dE* . (28) 

Taking into account Eq. (|2~TJ) . this last equation can be rewritten as follows: 

Pi = (P°)h + (a 2 /167rG) J [dj&ej + d 3 e t ) - 2d t d k e^ dY>\ , (29) 

where {P°)h is the energy corresponding to the first term of the r.h.s. of Eq. (|27|) . whose 
Fourier transform, /^(k), is given by Eq. fl25|) . As a result of the conditions satisfied by the 
quantities e»j (see the Appendix), it is easily proved that the term (P°)h vanishes. Hence, 

P T ° = (a 2 /lQnG) fidget - didft)dZ\ , (30) 

and, finally, the Gauss theorem allow us to write the surface integral in the last equation 
as a vanishing volume integral whatever the functions e^ may be and, in particular, for the 
functions leading to Eq. (|26l) . Hence, we have proved that the energy associated to any 
distribution of gravitational waves (propagating in a flat universe) vanishes; namely, the 
equation 

P T ° = 0. (31) 

is satisfied. 

Therefore, according to our protocol, we conclude that the energy P° due to tensor 
perturbations of a flat FRW vanishes. Notice that this conclusion does not depend either 
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on the spectrum P T {k) of the tensor perturbations (see the Appendix), or on the statistical 
character of the distribution of these perturbations. 

Notice again that similarly to what has been explained to the end of Sec. IHI A[ the 
resulting vanishing energy does neither depend on the chosen value of the to parameter, nor 
on the choice of S 3 . 

Since the energies due to scalar and tensor perturbations add, realistic perturbed universes 
including inflationary scalar modes have in all an infinite energy and so are not creatable. 
The presence of arbitrary tensor perturbation (zero energy) is irrelevant. Coming back to the 
end of Sec. IIII Al we see now that the energy of the scalar perturbed closed FRW universes is 
infinite irrespectively of, not only the rotation, traslation and dilatation groups, but actually 
on any infinitesimal coordinate transformation. 

IV. CONCLUSIONS AND DISCUSSION 

Our main conclusion is that perturbed flat FRW universes, including arbitrary tensor 
perturbations, and the adiabatic Gaussian scalar ones generated during standard inflation, 
have an infinite energy which is due to the scalar perturbations (see Sees. IIII Al and IIII Bj) . 



Since the total energy does not vanish, perturbed flat universes are not creatable, at least, 
in the framework of the standard inflationary paradigm, which appears to be compatible 
with most current observations. This conclusion implies that, among the perturbed FRW 
universes undergoing ordinary inflation, only the closed ones are creatable and, consequently, 
the slightly inhomogeneous universe where we live should be closed. 

It is generally believed that, in classical terms, there is no any way to decide if our universe 



is flat or closed, at least, if the curvature is small enough. An exception can be found in 15], 
where it is claimed that the spiral geodesic effect could be used to decide, observationally, 
whether we live either in a flat or a closed perturbed FRW universe. Another different 
method to distinguish between flat and closed perturbed universes arises from this paper. 
It is not directly based on observations. In our case, the creatable character of closed, and 
flat models compatible with observations, is studied according to the protocol described in 
Sec. HI The closed universes are creatable whatever the linear perturbations may be. This 
is a robust conclusion which privileges closed models against the flat ones. In the K = 
case, we have studied the most accepted model based on standard inflation. It does not 
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appear to be creatable; nevertheless, this conclusion is not valid whatever the perturbations 
may be. On the contrary, it is based on some assumptions and, consequently, it must be 
revised if some of such assumptions are modified in future. Particular attention deserve 
our hypothesis about: (i) the statistical isotropy of the universe and, (ii) the adiabatic and 
Gaussian character of any scalar perturbation in the post-inflationary era for k min < k < oo. 
Statistical isotropy has been recently questioned 19j [20|, and isocurvarture and (or) non 
Gaussian perturbations (based on cosmic strings, pre-inflationary processes and, so on) are 
not completely forbidden for k m i n < k < oo. Perturbed flat universes violating condition 
(i), and (or) condition (ii), and (or) any other possible condition, require particular studies 
if one wants to probe its creatable character. 



The subject of the creatable character of the 



considered in two papers in the first eighties 



6] 



non perturbed) FRW universes 
(there is also related work in [8] 



3 J 

,8 



ias been 
, In 



|6fl [7(1 , the authors discussed the possibility that the Universe could have arisen by quantum 
tunneling from 'nothing'. In [7j, a cosmological model is proposed in which the Universe 
is created by quantum tunneling from 'nothing' into a particular closed FRW model. Fur- 
thermore, in 6|, the authors find that within the context of FRW models, only the spatially 
closed and the flat de Sitter universes can originate in this manner, because they find that 
a finite tunneling amplitude exists only from initial spaces with finite three-volume (on the 
Euclidean section). 

In the absence of perturbations, the results in the present paper essentially agree with 
those of [6(-(9|, since we find with many other authors (including paper l|) that the closed 
and flat FRW models have vanishing energy and momenta and so, according to our termi- 
nology, this kind of universes would be creatable. Nevertheless, the method and the scopes 



ij-fl. 



of our work are very different from those of papers 

First of all, in these papers, the authors considered precise quantum mechanisms to orig- 
inate our classical universe from a quantum one; however, we assume that the produced 
classical universes must have vanishing energy and momenta and, then, we apply a definite 
protocol to decide whether a given classical universe (in the present case, some perturbed 
FRW models) has or not vanishing energy and momenta. Since we expect that any rea- 
sonable quantum process could not produce a universe with non vanishing energy, we have 
called these universes with vanishing energy and momenta creatable universes. In the above 
quoted papers, the authors assume that energy and momenta can only be properly defined 
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in asymptotic Minkowskian universes and, consequently, they could not follow our line of 
research. Nevertheless, in the present paper, we have been able to define the energy and mo- 
menta of some non asymptotic Minkowskian universes in a consistent and unambiguous way 
(the basic idea supporting our procedures). Asymptotically, perturbed and non perturbed 
FRW universes appear to be conformally flat, and this fact has allowed a definiton of energy 
and momenta which is a generalization of that used in the case of asymptotic Minkowskian 
space-times. In both cases, we are constrained to calculate the energy and momenta in ap- 
propriate cordinate systems making explicit a certain form of the spatial metric and, then, 
showing how cooordinate transformations preserving this form do not alter the energy and 
momenta. 

Furthermore, whereas in papers |G|] — 19|] only exact FRW universes were studied, we have 
considered realistic perturbed FRW universes (a related study considering departures from 
FRW geometries was proposed in [6|). In this way, we have found that perturbed closed FRW 
universes are creatable, but perturbed flat FRW universes, in the framework of standard 
inflation, are not. 

Notice that, attending the different criteria put forward to define a universe as creatable 
(to have a finite tunneling amplitude, in the case of the quoted authors, or to have vanishing 
energy and momenta, in our case), it is not obvious that both definitions must lead to the 
same conclusions. Comparison is possible for unperturbed FRW models (the quoted authors 
have not studied the perturbed ones). There is full agreement in the closed case, but for 
non perturbed flat models, there is some discrepancy, since the quoted authors find them 
creatable only in the de Sitter subcase, whereas we find all them creatable. Nevertheless, 
this discrepancy has perhaps a non significant meaning, since when we have considered a 
realistic perturbed flat FRW model, we have found that this model is not creatable. This 
result shows that the exact flat FRW universe is an unstable one to our concerning. In the 
case of perturbed FRW universes, the comparison with the results of these authors is not 
possible since, as we have said, they have not considered the case. 
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APPENDIX: SCALAR AND TENSOR INFLATIONARY PERTURBATIONS IN 
A FLAT BACKGROUND 

Our protocol begins with the use of Gauss coordinates, in which, the metric has the form: 

ds 2 = -dt 2 + dl 2 , dl 2 = g ij dx i dx j , (A.l) 

this means that, in Gauss coordinates, the conditions 

(7oo = -1, 9oi = (A.2) 

are satisfied. In the case of a FRW universe with scalar perturbations, conditions (1A.2J) 
define the so-called synchronous gauge. In this gauge, a detailed study about the evolution 



of scalar perturbations in a flat FRW universe can be found in reference 13J. The authors 
of that paper (Ma and Bertschinger) studied scalar perturbations in a rather general FRW 
flat universe containing baryons, cold dark matter (CDM), neutrinos, and radiation. They 
also studied the evolution in the longitudinal gauge. In the synchronous gauge, Ma and 
Berstchinger expanded the most general scalar metric perturbation as it is done in Eq. (j7j). 
This expansion only involves two arbitrary functions: r](k,r) and /i(k, r). 
In the longitudinal gauge one can write: 

ds 2 = a 2 (r)[-(l + 2i/j)dT 2 + (1 + 2<f))5 ij dx i dx j ] , (A.3) 

where only two arbitrary functions are necessary to describe the most general scalar pertur- 
bation. Fourier expansions of ip and <fi involve the coefficients ip(k, r) and 0(k, r). 

There are other coefficients appearing in the expansions of physical quantities involved 
in the energy momentum tensor, e. g., those corresponding to the density contrasts of the 
different energy components: <5 7 (k, r) for radiation, 5 c (k, r) for CDM and so on. All these 



coefficients are coupled in a complicate system of equations (see |13|). 

Fortunately, in order to evaluate the integrals giving the total energy and momenta of 
the perturbed universe, only functions f](k, r) and /i(k, r) (related to the metric) could be 
actually necessary; this fact facilitates our calculations. Moreover, only functions 7/(k, tq) 
and h(k, r ) are required, r being an arbitrary time. We can say that our problem is identical 
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to that solved in paper [13j, where initial conditions to solve the equations governing the 
evolution of the perturbations were calculated at a fixed time. This time was chosen to 
be in the time interval limited by electron positron annihilation and the time at which 
light massive neutrinos become non relativistic. The same choice is appropriate for us. Ma 
and Bertschinger solved the evolution equations in the mentioned time interval for linear 
adiabatic perturbations larger than the horizon (kr <<< 1). In the synchronous gauge, 
these authors found: 

V = 2C+ 5 + 4i ^ Ck 2 r 2 , h = Ck 2 r 2 , (A.4) 

' 6(15 + 4R U ) ' ' v ; 

where C = C(k), R u = p u /(p^ + f>u), and p v and p 7 are the background energy densities of 
neutrinos and photons, respectively. There is also a time decaying term whose explicit form 
is given in Eq. (T2U|) . 

The same study was also performed in the longitudinal gauge. Indeed, Ma and 
Bertschinger fixed their initial conditions by assuming that, in the longitudinal gauge, under 
the assumption of statistical isotropy, the power spectrum of the ip potential of Eq. (IA.3J) is 



P{if>) oc k~ 3 , (A.5) 

(see below for comments about this choice). How can we obtain the corresponding initial 
conditions in the synchronous gauge? The answer is easily obtained from Eqs. (18) in 13) . 
One of these equations reads as follows: 

^(k, r) = JL (h(k, r) + 677(k, r) + r^^k, r) + 6i)(k, r)]) , (A.6) 

where each dot stands for a derivative with respect to the conformal time. 
Taking into account Eqs. ()A.4l) - (lA.5j) and this last equation, one easily gets 



C(k) oc k~ 3/2 , (A.7) 

which is valid in the synchronous gauge. That could be our basic condition in order to 
compute the integrals giving the total energy and momenta of the universe. It is also the 
basic assumption leading to the initial conditions used by Ma and Berstchinger to solve the 
evolution equations in the synchronous gauge. The resulting numerical solution gave a very 
good description of both the power spectrum, P(k), of the energy density perturbations and 
the angular power spectrum (Ce coefficients) of the CMB. 
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Now, a question arises: why the ^-spectrum defined in Eq. ( 1A.5|) is appropriate? Let us 
try to answer this. 

After evolution, large enough cosmological inhomogeneities reenter the horizon in the 
matter dominated era and, afterward, it is well known that the potential if) satisfies the 
equation (see [13J |) 

Aifj oc 5 . (A.8) 

Since 5 is the total energy density contrast, this last equation indicates that function ip plays 
the role of the peculiar Newtonian gravitational potential. Finally, Eqs. (1A.5|) and (1A.8J) 
lead to 

P(k) = (\5 k \ 2 )^k, (A.9) 

which means that the spectrum of the energy density perturbations is a Harrison-Zel'dovich 
(HZ) one. This result justifies the use of the -^-spectrum defined in Eq. (IA.5I) . The above 
HZ spectrum is only valid at times close enough to horizon crossing, but afterward, as 
the inhomogeneities evolve inside the horizon, microphysics becomes important and this 
spectrum evolves toward a new one of the form P(k) = k/T(k), where T(k) is the so-called 
transfer function. 

Inflationary predictions are compatible with a HZ spectrum P(k) oc k, as well as with an 
spectrum of the form P(k) oc k ns having its spectral index n s close to unity. Accordingly, 
the analysis of the data obtained by the WMAP mission during three observation years 
leads to the inequality 0.942 < n s < 0.974. Moreover, if other observational data (galaxy 
surveys, other CMB observations, and so on) are taken into account, the resulting inequality 



appears to be 0.932 < n s < 0.962 (see [l6[ and [211]). From Eq. ( 1A.8I) one easily proves that 
the condition C(k) oc &( ns ~ 4 )/ 2 leads to a final spectrum P(k) = k ns . On account of these 
considerations, our calculations of the total energy and momenta of the universe is based on 
the relation: 

C{k) oc k {ns - 4)/2 , n s <l, n s ~l, (A.10) 

which coincides with Eq. ( 1A.7I) for n s = 1. 

We must emphasize that Eqs. ( 1A.4I) and ( 1A.10I) are only valid for adiabatic scalar pertur- 
bations evolving outside the horizon; hence, these relations only hold for k < k max , where 
kmax — 27r/L and L Q = H~ 1 (tq); here, L is the horizon size at the conformal time To- 
Moreover, if the adiabatic perturbations are inflationary, another inequality, k > k min , must 
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be also satisfied, where k min = 2ir/Lio and Ljo is the size, at time To, of a region compa- 
rable to the effective horizon at the beginning of inflation; that it to say, Ljo is the typical 
size of the huge inflationary bubbles. For k > k max the perturbations evolve inside the 
horizon (where microphysics is important) and, consequently, the spectra of super-horizon 
perturbations must be modified by means of transfer functions. If these perturbations are 
inflationary, they are initially Gaussian and afterward, during the radiation dominated era 



(in particular, in the period considered in paper [l3| and also in Sec. lIH A]) , they keep Gaus- 
sian because nonlinear processes leading to deviations from Gaussianity had not developed 
yet. Moreover, in the interval k m i n < k < oo, we assume that the perturbations produced 
during inflation are absolutely dominant against possible residual pre-inflationary fluctua- 
tions. The main reason is that, at the end of inflation, inflationary supercooling had made 
the pre-inflationary radiation density negligible against the total energy of the inflationary 
field and, consequently, after reheating, the mean radiation energy is fully dominated by 
the energy coming from the mean inflationary field and, evidently, the resulting adiabatic 
perturbations are associated to the fluctuations of this dominant field, with negligible contri- 
butions from pre-inflationary supercooled sources. Finally, for k < k m i n , the perturbations 
are so long that they will have a pre-inflationary origin without any inflationary contribu- 
tion. Thus, though these pre-inflationary perturbations can be expected to be small, their 
contribution to the integral in Eq. ffT2l) could be significant as k tends to zero and r tends 
to infinity. Thus, as a precaution, this interval has been also considered along the paper. 
Of course, it has been taken into account that these pre-inflationary perturbations could be 
non Gaussian. All these ideas are carefully taken into account in Sees. IIII Al and IIVI 

In general, inflation produces both scalar and tensor perturbations of the background 
universe. Some general considerations about tensor perturbations are now worthwhile. The 
tensor metric perturbations of a flat universe can be written in the form: 

^.(x,r) = /"d 3 fce ik - x /ij(k,r) = /"d 3 fce jkx i/ T (k,r)e y (k) , (A.ll) 

where functions e^- satisfy the following equations: 

€ij = €ji, €a = U, €ijKi = U , (A.izJ 

which ensure that quantities h^ (x, r) are symmetric, traceless, and divergenceless, as it must 
be in the case of metric perturbations describing gravitational waves. It is noticeable that 
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functions e^ only depend on the unit vector k and, consequently, any dependence on the 
scale (on k) of the tensor metric perturbation is involved in the coefficient H T (k,r). This 
scale dependence is usually fixed by defining a new power spectrum ( |22| ) 

P T (k) = k 3 (\H T (k)\ 2 ) oc k n T , (A.13) 

where n T is the tensor spectral index. With this spectrum, quantities ^(k) are proportional 
to /;;( n T~ 3 )/ 2 . Taking into account previous formulas and assumptions, we could calculate the 
total energy and momenta of the universe for different n T values. The explicit form of 
functions e^ is not necessary. Nevertheless, an explicit representation of these quantities 



can be easily found from the definitions given in [231 ] . What can we say about the spectral 
index n T ? 

The spectru m g enerated by most inflationary potentials has an spectral index n T ~ 1 — n s 



(see reference |24j); hence, from previous comments about n s -values and WMAP data, it 
follows that the power spectrum of these inflationary backgrounds of gravitational waves is 
very flat (small but non vanishing n T value). In the presence of these gravitational waves, 
the spectral index n s deviates from unity; by this reason and with the essential aim of 
allowing the existence of an inflationary background of gravitational waves (which is studied 
in Sec. IIIIBI) . the n s values used in Sec. IIII Al have been assumed to be slightly smaller than 
unity. 
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